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Abstract 

We derive Cayley's type conditions for periodical trajectories for the billiard 
within an ellipsoid in the Lobachevsky space. It appears that these new conditions 
are of the same form as those obtained before for the Euclidean case. We explain 
this coincidence by using theory of geodesically equivalent metrics and show that 
Lobachevsky and Euclidean elliptic billiards can be naturally considered as a part 
of a hierarchy of integrable elliptical billiards. 

MSC2000: 70H06, 53D25; PACS No 45.05,02.30H 

Subj. Class. Geometry of integrable systems 

Keywords: integrable billiards, Poncelet's theorem, Cayley's condition, spectral 
curve, geodesic hierarchy, separable perturbation 

1. Introduction 

We start with the following well-known integrable mechanical system: motion of 
a free particle within an ellipsoid in the Euclidean space of any dimension d. On the 
boundary, the particle obeys the billiard law. Integrability of the system is related 
to classical geometrical properties of elliptical billiards: the Chasles, Poncelet and 
Cayley theorems. According to the Chasles theorem [1] every line in this space is 
tangent to d — 1 quadrics confocal to the outer ellipsoid. Even more, all segments of 
the particle's trajectory are tangent to the same d—1 quadrics [26]. The Poncelet 
theorem [28, 22, 13] put some light on closed billiard trajectories: there exists a 
closed trajectory with d—1 given confocal caustics if and only if infinitely many 
such trajectories exist, and all of them have the same period. Since the periodicity 
of a billiard trajectory depends only on its caustic surfaces, it is a natural question 
to find an analytical connection between them and corresponding period. 

The Poncelet theorem, as one of the highlights of the XIX century projective 
geometry, attracted the attention of Arthur Cayley for several years (see [7-12]). 
In [8], Cayley found the analytical condition for caustic conies in the Euclidean 
plane case. The classical and algebro-geometric proofs of Cayley's theorem can 
be found in Lebesgue's book [28] and Griffiths and Harris paper [23], respectively. 
The generalisation is established by Dragovic and Radnovic for any d [18, 19]. This 
generalisation was done by use of the Veselov-Moser discrete quadratic L — A pair 
for the classical Heisenberg magnetic model [32]. 
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The integrability of elliptical billiard systems in the Lobachevsky space was 
proved by Veselov in [38]. There, Veselov used discrete linear L — A pair, which is 
quite different from the one used in the Euclidean case. 

The starting point of this paper is derivation of Cayley's type conditions for the 
Lobachevsky billiard and our observation that these new conditions coincide with 
those obtained in [18, 19] for the Euclidean case (Section 3). 

We found a natural way to explain this coincidence and it is related to the 
recently developed integrability approach in the theory of geodesically equivalent 
metrics [29, 35]. Both Lobachevsky and Euclidean elliptic billiards can be naturally 
considered as members of a hierarchy of integrable elliptical billiards (Section 4). In 
the conclusion of this Section, we present some properties of the Laurent polynomial 
integrable potential perturbations of those separable systems, continuing the study 
of such systems which started with [14], see also [15, 24, 17, 16]. 

2. Basic notions on billiard systems 



Let (Q, (?) be a d-dimensional Riemannian manifold and let D C Q be a domain 
with a smooth boundary F. Let tt : T*Q ^ Q be a natural projection and let 
be the contravariant metric on the cotangent bundle, in coordinates 



Consider the reflection mapping 

r : 7r~"^r — > 7r~"^r, p_ i— > p+, 

which associates the covector p+ e T*Q, a; e F to a covector p_ e T*Q such that 
the following conditions hold: 



(1) 



\p+\ = \p-U 
p+-p-±r. 



A billiard in D is a dynamical system with the phase space M = T*D whose 
trajectories are geodesies given by the Hamiltonian equations 

dH . dH , \ , 

(2) "^"a^' ^(P'^) = 2^=^ 

reflected at points a; e F according to the billiard law: r(p_) = p+. Here p- and 
p+ denote the momenta before and after the reflection. If some potential force field 
V{x) is added than the system is described with the same reflection law (1) and 
Hamiltonian equations (2) with the Hamiltonian H{p,x) = ^g~^{p,p) + V{x). 

A function / : T*(5 ^ M is an integral of the billiard system if it commutes with 
the Hamiltonian {{f,H} = 0) and does not change under the reflection {f{x,p) = 
f{x,r{p)), X e F). The billiard is completely integrable in the sense of Birkhoff if 
it has d integrals polynomial in the momenta, which are in involution, and almost 
everywhere independent (see [26]). 

The classical integrable examples, with smooth boundary, are billiards inside 
ellipsoids on the Euclidean and hyperbolic spaces and spheres, with integrals qua- 
dratic in the velocities [26]. These systems can be also considered as discrete in- 
tegrable systems [37, 38] . The explicit integrations in terms of theta-functions are 
performed by Veselov, Moser and Fedorov (see [37, 32, 38, 20]). 
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3. PONCELET THEOREM AND CAYLEY'S CONDITION 
FOR THE BILLIARD IN THE LOBACHEVSKY SPACE 

Vcsclov proved the integrability of the bilhard system within an ellipsoid in 
the Lobachevsky space in [38]. He showed that its motion corresponds to certain 
translations of the Jacobi variety of some hyperelliptic curve and gave explicit 
formulae of the motion in terms of thcta-functions. The aim of this section is to 
find an analogue of Poncelet's and Cayley's theorem [8] for the billiard motion 
within an ellipsoid in the Lobachevsky space. 

3.1. Integration of the billiard motion in the Lobachevsky space. Pon- 
celet theorem. For a brief account of Vesclov's results on the billiard in the 
Lobachevsky space [38], let us consider the {d + l)-dimensional Minkowski space 
V = R**'-*^ with the symmetric bilinear form: 

V) = -^oVo + ^iVi H \- UVd- 

One sheet of the hyperboloid ^) = — 1 with the induced metric is a model of the 
d-dimensional Lobachevsky space H''. An ellipsoid F in this space is determined by 
the equation 

(3) r= l^eH'^, -^ + ^ + --- + ^ = ol, 

I ao ai aa J 

with ao > ai > 02 > • • • > fld > 0. 

All segments of the billiard trajectory within this ellipsoid are tangent io d — 
1 confocal quadric surfaces (including multiplicity), fixed for a given trajectory 
(Theorem 3 in [38]). Denote hy Hi, i = 1, . . . ,d — 1 the numbers such that the 
equations of these caustics are: 

11 1 

(4) ^ + + + =0, (l<i<d-l). 

ao - in ai - /ii ad - m 

Then the points of reflection from the boundary F correspond to the shift Dk+i = 
Dk + Q- — Q+ on the Jacoby variety of the spectral curve C 

(5) C : (/i-ao)...(/U-ad) =c- A^(/i-/ii)...(/i-/id_i), 

where c is a constant, and Q+, Q- are the points on the curve C over n = 0. (See 
Theorem 2 of [38] . The curve C is the spectral curve of the L — A pair considered 
there.) 

Let us note that Veselov considered only the case of the regular (hyperelliptic) 
curve C [38]. However, his consideration holds for the singular case, too. 

Suppose a periodical billiard trajectory inside the ellipsoid F in the Lobachevsky 
space is given. All trajectories with the same caustics have the same spectral curve. 
If the period of the given trajectory is n, then n((5+ — Q-) = on Jac{C), and 
vice-versa. Thus, all these trajectories close after n bounces. Therefore, Poncelet's 
-type theorem for the billiard in the Lobachvsky space is derived from Veselov's 
results: 
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Proposition 1. Suppose a periodical billiard trajectory inside an ellipsoid in the 
Lobachevsky space is given. Then any billiard trajectory which shares the same 
caustic quadrics is also periodical, with the same period. 

3.2. Cayley's conditions — regular spectral curve. Assume that all con- 
stants fli, • • • , Old, A*i, • • • , A*d-i are mutually different. Then the spectral curve 
C is hypcrcUiptic. Cases when some of them coincide are discussed in the next 
subsection. 

To establish an analytical condition on a trajectory to be periodic with period 
n, wc need to find out when the divisors nQ+ and nQ- on the spectral curve are 

equivalent. 

Lemma 1. Let the curve C be given by 



(6) 



y 



{X-Xi)...{x- X2g+2), 



with all Xi mutually different and not equal to 0, and Q+, Q- the two points on C 
over the point x = Q. Then nQ+ = uQ^ is equivalent to: 



(7) 



rank 



g+2 
's+3 



9+3 
S+4 



Bn+1 
Bn+2 



B. 



g+n 



B- 



2n-l 



< n — g and n> g, 



where y = \/{x — x{) ... (a; — X2g+2) = Bq + Bix + B2x'^ + . . . is the Taylor expan- 
sion around the point Q- . 

Proof. C is a hyperelliptic curve of genus g. The relation = nQ- means that 
there exists a meromorphic function on C with a pole of order n at the point (5+, a 
zero of the same order at Q- and neither other zeros nor poles. Denote by L(n(5+) 
the vector space of meromorphic functions on C with a unique pole of order 
at most n. Since is not a branching point on the curve, dimL(n(5+) = 1 for 
n < g, and dim i(nQ+) = n — g + 1, fov n > g. In the case n < g, the space 
L{nQ+) contains only constant fimctions, and the divisors nQ+ and nQ- can not 
be equivalent. If n > g + 1, we choose the following basis for L(nQ+): 



where 



fk = 



f) /l) • • • ) fn—gi 

y- Bp - Bix Bg+k-ix 

^g+k 



g+k-1 



Thus, nQ+ = nQ- if there is a fimction / G L{nQ^) with a zero of order n at 
Q_, i.e., if there exist constants ao, . . . , a„_g, not all equal to 0, such that: 

ao + ai/i(Q_) + ... Q;„-g/„-g(Q_) = 



ai/{(Q-) 



+ 



J'iQ-) = 



ai/l"-^)(Q_) + 



an 



-gftg'HQ-) = 0. 
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Existence of a non-trivial solution to this system of linear equations is equivalent 
to the condition (7). □ 

Introducing new coordinates x = ji, y = \/c\{ijl — fxi) . . . {h — l-id-i), the spectral 
curve (5) is transformed to: 

(8) y"^ = (x - ao) . ■ .{x - ad){x - in) ...{x- Hd-i), 

and wc obtain the following 

Theorem 1. The condition of a billiard trajectory inside the ellipsoid (3) in the 
d-dimensional Lobachevsky space, with non-degenerate caustics (4), to be periodic 
with period n> d is: 



rank 



Bn+l Bn ■ . . Bd+l 

Bn+2 Bn+l ■ ■ ■ Bd+2 



B'2n-^ B^n-?. ■ ■ ■ B. 



n+d-\ 



< n - rf+ 1, 



where \/ {x — ag) . . .{x — ad){x — ^i) . . . (x — A'd-i) = Bq-\-BiX-\-B2x'^ + There 

is no such trajectories with period less than d. 

3.3. Cases of singular spectral curve. When all ao, oi, . . . , a^, /ui, . . . , /id-i 
are mutually different, then the curve (5) has no singularities in the afEne part. 
However, singularities appear in the following three cases and their combinations: 

(i) ai = iij for some i,j. The spectral ciuvc (5) decomposes into a rational and a 
hyperelliptic curve. Geometrically, this means that the caustic corresponding to jii 
degenerates into hyper-plane Xi = 0. The billiard trajectory can be asymptotically 
tending to that hypcr-planc (and therefore can not be periodic), or completely 
placed in this hyper-plane. Therefore, the closed trajectories appear when they are 
placed in a coordinate hyper-plane. Such motion can be discussed like in the case 
of dimension d — 1. 

(ii) Qi = Qj for some i ^ j. The ellipsoid (3) is symmetric. 

(iii) Hi = fjtj for some i ^ j. The billiard trajectory is placed on the correspond- 
ing confocal quadric hyper-surface.* 

In the cases (ii) and (iii) the spectral curve C is a hyperelliptic curve with sin- 
gularities. In spite of their different geometrical nature, they both need the same 
analysis of the condition nQ+ = nQ- for the singular curve (5). 

Lemma 2. Let the curve C be given by y^ = {x — xi) ... (a; — 0:23+2), with all 

Xi different from 0, and Q^, the two points on C over the point x = Q. Then 
nQ+ = nQ- is equivalent to (7) where y = ^ {x — xi) . . .{x — X2g+2) = Bq + Bix + 
B2x'^ + . . . is the Taylor expansion around the point Q_ . 

Proof. Suppose that, among a;i, . . . , X2g+2, only X2g+i and X2g+2 have same values. 
Then (.T2g+i,0) is an ordinary double point on C. The normalisation of the curve 
C is the pair (C, tt), where C is the curve given by: 

C -.y^ = {x-xi)...{x- X2g), 



*We learned about its geometric significanse from Prof. Yurij Fedorov 
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and TT : (5 — > C is the projection: 

{x,y) ^{x = x, y = {x- X2g+i)y). 

The genus of (7 is (/ — 1. The relation nQ+ = nQ- is equivalent to existence of 
a meromorphic function / on (7, / G L{nQ+), with a zero of order n at Q-, and 
f{A) = f{B), where Q+, Q- are the two points over x = 0, and A,B are over 

X = X2g+1- 

For n< g—1, dim L{nQ+) = 1, and this space contains only constant functions. 
For n> g, we can choose the following basis for L{nQ+): 

l,/o,/i°7r, ...,/„_g0 7r. 

fk are as in Lemma 1 for A; > 0, and 

, _y-Bo-Bix Bg_ixB-^ 



where y = ^{x — xi) . . . {x — X2g) = Bo + Bix + B2SP' + . . . is the Taylor expansion 
around the point Q-- 

Since /o is the only element of the basis with different values in the points A and 

we obtain that nQ^ = nQ^ is equivalent to (7). 

Cases when C has more singularities, or singular points of higher order, can be 
discussed in the similar manner. □ 

Immediate consequence of Lemma 2 is that Theorem 1 can be applied not only 
for the case of the non-singular spectral curve, but in the cases (ii) and (iii) too. 

Thrcrcforc, the following interesting property holds. 

Theorem 2. If the billiard trajectory within an ellipsoid T in d- dimensional Loha- 
chevsky space is periodic with period n < d, then it is placed in one of the n- 
dim,ensional planes of symmetry of the ellipsoid. 

This property can be seen easily for d = 3. 

Example 1. Consider the billiard motion in an ellipsoid in the 3-dimensional 
space, with /zi = /Lt2, when the segments of the trajectory are placed on generatri- 
ces of the corresponding quadric surface confocal to the ellipsoid. If there existed 
a periodic trajectory with period n = d = 2>, the three bounces would have been 
complanar, and the intersection of that plane and the quadric would have consisted 
of three lines, which is impossible. It is obvious that any periodic trajectory with 
period n = 2 is placed along one of the axes of the ellipsoid. So, there is no periodic 
trajectories contained in a confocal quadric surface, with period less or equal to 3. 

4. Hierarchy of integrable elliptical billiards 

4.1. The Beltrami— Klein model of the Lobachevsky space. Note first that 

Caylcy's type conditions for the Lobachevsky billiard from Section 3 are of the same 
form as those obtained in [18, 19] for the Euclidean case, although the L — A pairs 
used there are quite different. There is a natural way to explain this coincidence. 
We use the Beltr ami-Klein model of the Lobachevsky space H''. 
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The coordinate transformation 

^1 u 

yi = -r^---^yd = -r 
?o to 

maps the Lobachevsky space, modeled as a pseudosphere of the Minkowski space, 
to the Beltrami-Klein model within the unit sphere in R'' [38] . 

Now, after appropriate linear changing of coordinates xi = aiyi, . . . ,Xd = otdUd 
we can obtain the Beltrami-Klein model inside the ellipsoid A: 

'^={-=(^' "'^^^'•vyh-^vr']- 

such that the ellipsoid (3) in new coordinates is confocal to A. Then its equation 
can be written in the form: 

C /y>2 ™2 

r = J a; e + ' ' ' ' 



' 6i — c 62 — c bd — c 

where < c < i = 1, . . . ,d. 

The hyperbolic metric within A is given by (for example, see [34]): 

_ - , 4^ f (^^1 , , dXd\ , f -El^fi , .Xddxd 

- - - ^ -'H&i bdJ^K b, bd , 




bi ' b2 ' ' bd 



1 , 2 ^ 



The metric dg"^ can be written in the matrix form as dg"^ = {Ildx,dx), where 
dx = {dxi, . . . , dxd), 

n=d^(^^"+^"^^^"^)' 

B = diag(6i, . . . , bd), and (•, •) is the Euclidean scalar product. 

The hyperbolic metric has the same geodesies, considered as unparemetrized 
curves, as the Euclidean metric dg"^ = dx^ + dx^ H + dx^^. 

Suppose that 61 > 62 > • • • > 6^. The standard elliptic coordinates Ai, . . . , Ad in 
M'' (61 > Ai > 62 > A2 > • • • > A<j_i > bd > Ad) are defined as solutions of the 
equation: 

22 2 

The direct verification shows that the metric; d.<f , as well as the Euclidean metric 
dg^, is orthogonally separable in the elliptic coordinates and geodesic flows can be 
integrated by the theorem of Stackel. This means that hypersurfaces Aj = const of 
the coordinate system Ai, . . . , A^ are orthogonal to each other and the corresponding 
Hamilton- Jacobi equations for the Hamiltonian of the geodesic flows have complete 

solutions of the form 5(Ai, . . . , Ad, ci, . . . , q) = 5i(Ai, ci) H |-<S„(Ad, Cd) (see [2, 

3] and references therein). 
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Consider the billiard in the domain D bounded by the ellipsoid T. In elliptic 
coordinates, the boundary of the ellipsoid is given by the equation Xd = c and the 
reflection map, both for the Euclidean and Lobachevsky metrics, is given by 

(Al,A2,...,Ad-l,C,PAi,PA2,---,PAd_i,PAj ^-> 

(Ai, A2, . . . , Ad-1, c,pAi,PA2, • • • ,PAd_i, -PaJ, 

where (A, px) are canonical coordinates in r*M''. Such a simple form of the reflection 
map is due to the fact that F and A are confocal in the coordinates xi,. . . ,Xd- 
Therefore we have 

Lemma 3. The billiards inside the ellipsoid T in the Euclidean and the Loba- 
chevsky space, modeled within the ellipsoid A, have the same trajectories up to 
reparametrization. 

Lemma 3 provides the explanation for the coincidence of the Cayley's conditions 
obtained in the previous section and papers [18, 19]. The above observation allows 
us to approach to the problem of the integrability of elliptical billiards in a new 
way, using theory of gcodesically equivalent metrics. 

4.2. Geodesically equivalent metrics. Let g and g be Riemannian metrics on 

d dimensional manifold Q. The metrics g and g are called geodesically equivalent 
if they have the same geodesies considered as unparemetrized curves. This is a 
classical subject studied by Beltrami, Dini, Levi-Civita, etc. in 19th century. Re- 
cently, the new global uncrstanding of the theory is developed in the framework of 
integrable systems (see [35, 36, 29, 5] and references therein). 

Having the metrics g and g, deflne the (l,l)-tensor field L = L{g,g) by 

Consider functions 

(10) Mp, x) = g-\Si p,p) = J2 {Si)W''PiPk' 

where (1,1) tensors Sfe are given by the formula: 

{det(L+aId)}(L+aId)-i = Sd_i a'^''^ + Sd_2 a**"^ + • • • + So, a e M 

If the metrics g and g are geodesically equivalent then functions Ji (p, x) are in 
involution with respect to the canonical symplectic structure on T*Q. Moreover, 
if the eigenvalues of L arc all different at one point of Q, then they arc different 
almost everywhere and the geodesic flows of g and g are completely integrable. The 
complete set of involutive integrals for the first fiow is Jo, Ji, . . . , Jd-i (see [35, 29, 
36]). In this case we say that g and g are strictly nonproportional. 

The pair g, g of geodesically equivalent Riemannian metrics produces the family 
of gcodesically equivalent Riemannian metrics gk, gk, k G Z given by the formulas 
[35, 36]: 

(11) 5fe(C,??)=5(L'C,??), 5fe(C,??)=5(L'=C,??), keZ. 
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The integrals of the geodesic flows of metric gk and gk+2 canonically given by (10) 
coincide [35]. Following [5] we call (11) Topalov-Sinjukov hierarchy of Riemannian 
metrics on Q. 

The nice geometrical interpretation of geodesical equivalence is done by Bolsinov 
and Matveev [5] . They proved that g and g are geodesically equivalent if and only 
if L is a Benenti tensor field for the metric g [5, 3]. This implies that if g and g are 
strictly nonproportional than all metrics gk, gk are orthogonally separable in the 
same coordinates and geodesic flows can be integrated by the theorem of Stackel 
[5]. 

4.3. Hierarchy of integrable elliptical billiards. Now we shall apply the 
general construction described in the previous section to the Euclidean dg^ and 
Lobachevsky metrics dg^ inside the ellipsoid A. Note that this natural geodesical 
equivalence is a slight modification of the geodesical equivalence studied by Topalov 
in [36]. Taking b — B, d — ^/^x, from lemma 7 of [36] we are getting that (1,1) 
tensor field L has the following matrix form 

L = L{dg^,df) = (det n)^ H"^ = B -x x. 

Therefore we have the Topalov-Sinjukov hierarchy of strictly nonproportional 
Riemannian metrics within A given by: 

dgl = {{B-xi^x)'^dx,dx), 
dgl = (n(B -X (g) x)''dx, dx). 

The corresponding geodesic flows are completely integrable. The integrals of the 
geodesic flow of the metrics dgl and dg1_^2 =^^6- 

(13) J,'{p, x) = (Si(B -X x)-''p,p) 

where 

{det(B-a;(g)a; + aId)}(B-a;0a; + aId)~^ = 
^^"^^ = det B„ ((1 - (B-i x, x)) B-i + B-i B'^) = S^-i a"^"' + • • • + So, 



p = {p\,...,Pd) e T*'W^ is the canonical momentum in Euclidean coordinates 
x = {xi, . . . , Xd), Bq, = diag(6i + a, . . . ,bd + a) and a is a real parameter. For 
k = these functions are defined on the whole T*]!*^ and coincide with commuting 
functions given by Moser in [31]. 

Consider the billiards in the domain D bounded by the ellipsoid T with metrics 
(12). According to [5] all metrics (12) are orthogonally separable in elliptical co- 
ordinates. This implies that the reflection map is the same for all metrics and in 
elliptic coordinates has the form (9). Moreover, since integrals Jf (p, x) are diagonal 
in elliptic coordinates, we have that they arc not just integrals of the geodesic flows 
of metrics dgl and dgl_^_2, but also integrals of the corresponding billiard systems 
inside ellipsoid T (see [26], pages 133-134). 

Thus we get the following general statement: 
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Theorem 3. The billiard systems inside ellipsoid T with the Riemannian m,etrics 
dg^, dgl given by (12) are completely integrable for all k & Z. In particular, the 
elliptical billiards in the Euclidean and hyperbolic spaces are completely integrable. 

Since the reflection map r is the same for the whole hierarchy, applying Propo- 
sition 3 of [35], we get 

Corollary 1. The billiard systems for the given hierarchy have isomorphic Liou- 

ville foliations ofT*D/r. 

Remark 1. Matveev and Topalov [29] and Tabachnikov [34] proved that ellipsoid 

admits nontrivial geodesic equivalence between the standard metric and the metric 
1 / dxl dxl dx^^-^ 

6i-6,...6.+i.(fi + |i + ... + ;|±i) + + + 



lA- 



The Euclidean and the Lobachevsky metrics within ellipsoid A can be seen as limits 

of the given metrics as bd+i tends to zero. 

4.4. Integrable potential perturbations. We shall say that the potential V{x) 
is separable in the elliptic coordinates Ai, . . . , if the Hamilton-Jacobi equation 
for the Hamiltonian ^ (Pi + • • • ^Pd) + ^(■^) can be solved by separation of variables 
in elliptic coordinates. This definition, in a more geometrical fashion, can be found 
in [2, 3]. The potential of the elastic force is an example. 

The potential V{x) is separable in the elliptic coordinates on if and only if 
V{x) is a solution of the linear system of partial differential equations 

n . . d'^y f d d \ / , ^ dv\ ^ 



fe=i 



'dxk J 



for i ^ j (see [30, 2]). 

We shall denote the Hamiltonians of the geodesic flows of metrics dg^ and dg^ by 
Hq{p,x) and Hq{p,x) respectively. Suppose that V{x) is a solution of (15). Then 
from [5] follows that Hamiltonian systems with Hamiltonian functions 

(16) H'^ip, x) = H^ip, x) + V{x), H\p, x) = H^{p, x) + V{x) 

are completely integrable, and can be solved by separation of variables in elliptic 
coordinates for all k. There is a complete set of commuting integrals of the form 

I^{p, x) = J^{p, x)+fi{x), i = l,...,d, 

for each Hamiltonian H^{p^x) where J^{p,x) is given by (13). The functions /,;(.x) 
do not depend of k. They are solutions of the equations V fi{x) = SiVl^(a;), 
where V/ = {dif,..., dnf) and Si are given by (14). Similar statement holds for 
Hamiltonian systems with Hamiltonians H^(p,x). 

Consider the billiard systems with Hamiltonians (16) within the ellipsoid T. 
Prom the choice of Ii{p,x) we have that these functions do not change under the 
reflection. Thus we get the following corollary. 
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Corollary 2. Suppose that V{x) is a solution of (15). Then the billiard systems 
with Hamiltonians (16) within the ellipsoid T are com,pletely integrable. 

Let us consider the solution of equations (15) in the form of Laurent polynomials 
(17) V{x)= Piu-,iA4' ■■■^it, k-,k+GZ. 

k- <ii ,...,id<fc-|- 

Suppose that Laurent polynomial (17) is a solution of (15). Then coefficients 
Pii,...,id satisfy the following system of difference equations 

ih -bi)ikiiPn,...,i4 = 
^ ) {■ ■ \c ■ \ 

\ll + • • • + id)(i;Pii,...,i;s_i,i|s— 2,ife+i,...,i<j ~ *fePii,...,ii_i,ii— 2,ii+i,...,i<j)- 

Such potential perturbations are described for d = 2 in [14] (see also [15, 24]) 
and for d = 3 in [17]. In general, the linear space of Laurent polynomial solution 

of (15) has a basis of the form 

Vfe = Vk{xi, . . .,Xd), 
^^^^ Wl = ^Pi_,{x^,...,Xd), i = l,...,d, 

where Vfe and VI are polynomials of degree 2k, k> 0. The potentials Vfe and W^, 
in elliptic coordinates, correspond to the potentials 

with v{t) = Y!j=i ajf'''^'^^ and v{t) = l3j{t - ai)-^ respectively 
Example 2. As an example, we write down a few of the basis potentials (19): 



Vx{x) = x^j (Jacobi), V2(x) = ^i^i " 51 1 ' 

3 3 \ 3 

V3(x) = E ^^3 - 2 in f E + 



3 



1 I x"^- \ 

Wlix) = — (Rosochatius), >V^(a:) = — 1 + V 

Xi y bi bj J 



2 Av»2 ™2 \ ^2 ™2 

b, ' {b~b;yl ib,-b,r{b,-b,y 



Let V{x) = ^pOipVp{x) be some separable polynomial potential. Consider bil- 
liard systems with Hamiltonians (16). By the Maupertiues principle [1], for a given 
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value of total energies h, satisfying condition h > maxxeo ^(2^)1 the motions in the 
potential field V{x) inside T are reduced to geodesical motions with metrics 



(20) ih-Vix))dgl {h-V{x))dgl 

It is clear that the billiard systems within F with metrics (20) are integrable. 

Concluding remarks. Theorem 3 holds also if the boundary of the billiard is the 
union of the confocal quadrics T — UFc^ • • -UFc^, — {x £ R'', 7(0^) = 1}, or 
more generally, if the billiard is constrained to some of confocal quadrics. The same 
results can be formulated for billiards constrained on spheres by using geodesical 
equivalences established in [35, 29]. Then systems are orthogonally separable in the 
spherical elliptic coordinates. 

Polynomial potentials separable in elliptic coordinates on M.'^ and spheres S"^ are 
given by Bogoyavlenski [4] and Wojceichowski [39]. After Rosochatius's potential 
V{x) = Oi'i^i'^ (see [33] and Apendix in [31]), particular examples of rational 

potentials are found by Braden and Wojceichowski [6, 39]. Kalnins, Benenti and 
Miller described separable rational potentials in terms of certain recurrence relations 
between potentials of different degrees [25] . 

Recently, Fedorov integrated the elliptical billiard with the elastic potential in the 
Euclidean spase [21]. Dragovic found the connection between the Laurent potential 
perturbations of the elliptical billiards for rf = 2 and the Appell hypergeometric 
functions [16]. Also recently, the two-dimensional billiards with smooth boundary 
and additional irreducible integrals of the third and fourth degree, with a help of the 
integrable cases of Goryachev-Chaplygin and Kovalevskaya are given by Kozlova 
[27]. 
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